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It is proposed an integral formulation of classical Yang-Mills equations in the presence of sources, 
based on concepts in loop spaces and on a generalization of the non-abelian Stokes theorem for two- 
form connections. The formulation leads in a quite direct way to the construction of gauge invariant 
conserved quantities which are also independent of the parameterization of surfaces and volumes. 
Our results are important in understanding global properties of non-abelian gauge theories. 
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The aim of the present paper is to propose an inte- 
gral formulation of the classical equations of motion of 
non-abelian gauge theories. Our approach is based on a 
generalization of the non-abelian Stokes theorem for two- 
form connections, which allows to present the Yang-Mills 
equations as the equality of an ordered volume integral 
to an ordered surface integral on its border. The formu- 
lation leads in a quite simple way to the construction of 
gauge invariant conserved quantities which are indepen- 
dent of the parameterizations of volumes and surfaces. 
The most appropriate mathematical language to phrase 
our results is that of generalized loop spaces. There is 
a quite vast literature on integral and loop space for- 
mulations of gauge theories [1]. Our approach differs in 
many aspects of those formulations even though it shares 
some of the ideas and insights permeating them. We 
make however concrete progress in relation to those ap- 
proaches. The main statement of this paper is: 

Consider a Yang-Mills theory for a gauge group G, 
with gauge field A^, in the presence of matter currents 
J p , on a four dimensional space-time M . Let £1 be any 
tridimensional (topologically trivial) volume on M , and 
dfl be its border. We choose a reference point xr on d£l 
and scan Q with closed surfaces, based on xr, labelled by 

and we scan the closed surfaces with closed loops based 
on xr, labelled by t, and parametrized by a, as we de- 
scribe below. The classical dynamics of the gauge fields is 
governed by the following integral equations, on any such 
volume fl, 
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where Pi and P3 means surface and volume ordered inte- 
gration respectively, is the Hodge dual of the field ten- 
sor, i.e. Fpv = \ e^ p x F pX , with F M „ = d ih A v - d v A^ + 
i e [ A^ , A v \, e is the gauge coupling constant, a and (3 
are free parameters, and where we have used the notation 
X w = W^ 1 X W, with W being the Wilson line defined 
on a curve T, parameterized by a, through the equation 
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where x^ (fi = 0,1,2,3) are the coordinates on the four 
dimensional space-time M . The quantity V is defined on 
a surface S through the equation 
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with T (A, t) = ie / daW' 1 
and where 
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where Jp V \ is the Hodge dual of the current, i.e. J p — 
_ £ ni/p\ j vpX The Yang-Mills equations are recovered 
from (1) in the case where Q is taken to be an infinitesi- 
mal volume. Under appropriate boundary conditions the 
conserved charges are the eigenvalues of the operator 

Qs = P 2 e ie Ios dTda{aF " + ^^^ = P 3 eIs d<:dTVJV " 

(5) 

where S is the 3- dimensional spatial sub-manifold of M. 
Equivalently the charges are TrQg. 

In order to prove that (1) does correspond to and inte- 
gral formulation of the classical Yang-Mills dynamics, we 
shall start by describing the generalization of the non- 
abelian Stokes theorem as formulated in [2, 3], Con- 
sider a surface E scanned by a set of closed loops with 
common base point xr on the border <9E. The points 
on the loops are parameterized by a £ [0, 2tt] and each 
loop is labeled by a parameter r such that t = corre- 
sponds to the infinitesimal loop around xr, and r = 2ir 
to the border 9E. We then introduce, on each point 
of M, a rank two antisymmetric tensor B^ v taking val- 
ues on the Lie algebra Q of G, and construct a quan- 
tity V on the surface E through (3), but with T(A,t) 
replaced by T (B, A, r) = da W' 1 B^W^r^, 
and where the cr-integration is along the loop Y labeled 
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by r, and W is obtained from (2), by integrating it along 
r from the reference point xr to the point labeled by 
a, where is evaluated. By integrating (3), from the 
infinitesimal loop around xr to the border of E, we ob- 

tain V = Vr P 2 eJo Jo d ° d ^ , where P 2 

means surface ordering according to the parameteriza- 
tion of £ as described above, and Vr is an integration 
constant corresponding to the value of V on an infinites- 
imal surface around xr. If one changes £, keeping its 
border fixed, by making variations fa M perpendicular to 
E then V varies according to (see sec. 5.3 of [2], sec. 2.3 
of [3], or the appendix of [4]) 
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W- 1 [D X B^ + D^B„ X + D v B Xlx ] W^- ~ Sx x 

do UT 

dx K dx>* 
da 



£ da' [ BZ (a') - ieFZ (a') , B™ (a) ] ^ 

( d -^5x»(o)-6x P (o>) ^M)W-(r) 
\ dr dr J J 



where D^* = * 



i e [ An , * ] . The quantity V (r) ap- 
pearing on the r.h.s. of (6) is obtained by integrating 
(3) from the infinitesimal loop around xr to the the loop 
labelled by r on the scanning of E described above. Note 
that the two cr-integrations on the second term on the 
r.h.s. of (6) are performed on the same loop labelled by 
r. Consider now the case where the surface E is closed, 
and the border of E is contracted to xr. The expression 
(6) gives then the variation of V when we vary E keeping 
xr fixed. Therefore, if one starts with an infinitesimal 
closed surface T,r around xr one can blows it up until it 
becomes E. One can label all those closed surfaces using 
a parameter ( g [0,27r], such that ( — corresponds to 
T,r and ( = 2ir to E. The expression (6) can be seen as 
a differential equation on £ defining V on the surface E, 
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-JCV = 
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where JC corresponds to the r.h.s. of (6) with 6x^ re- 
placed by By integrating (7) from E^ to E, one 
obtains V evaluated on E, which is now an ordered vol- 
ume integral, over the volume 51 inside E, and the order- 
ing is determined by the scanning of 51 by closed surfaces 
as described above. But this result has of course to be 
the same as that obtained by integrating (3) when the 
surface is closed, namely <951. Therefore, we obtain the 
generalized non-abclian Stokes theorem for a two-form 
connection B^, parallel transported by a one-form con- 
nection A„ 
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where P3 means volume ordering according to the scan- 
ning described above, and Vr is the integration constant 
obtained when integrating (3) and (7). It corresponds in 
fact to the value of V at the reference point xr. Note 
that such theorem holds true on a space-time of any di- 
mension, and since the calculations leading to it make no 
mention to a metric tensor, it is valid on flat or curved 
space-time. The only restrictions appear when the topol- 
ogy of the space-time is non-trivial (existence of handles 
or holes for instance). 

Going back to (1) one notes that it can be ob- 
tained from (8) by replacing B^ by ie a Fp V + f3 F^ , 
and using the Yang-Mills equations, D v F vl1 — and 
D U F»» = 0, to replace (D X B^ + D„B vX + D v B Xfl ) in 
(6) by (— ie/3 J^ux), and so JC introduced in (7) is now 
given by JC — J Q 2 n dr V J V~ x , with J given in (4). 
Therefore, (1) is a direct consequence of the Yang-Mills 
equations and the Stokes theorem (8). Note that Vr in- 
troduced in (8), does not appear in (1) because it has 
to lie in the centre Z (G) of G to keep the gauge co- 
variance of (1) (see [5]). On the other hand the integral 
equation (1) implies the local Yang-Mills equations. In 
order to see that, consider the case where 51 is a infinites- 
imal volume of rectangular shape with lengths dx^ , dx v 
and dx x along three chosen Cartesian axis labelled by fi, 
v and A. We choose the reference point xr to be at a 
vertex of 51. By considering only the lowest order contri- 
butions, in the lengths of 51, to the integrals in (1), one 
observes that the surface and volume ordering become 
irrelevant. We have to pay attention only to the orien- 
tation of the derivatives of the coordinates w.r.t. the 
parameters er, r and £, determined by the scanning of 51 
described above. In addition, the contribution of a given 
face of 51 for the l.h.s. of (1) can be obtained by evalu- 
ating the integrand on any given point of the face since 
the differences will be of higher order. Consider the two 
faces parallel to the plane x^x v . The contribution to the 
l.h.s. of (1) of the face at xr is given by —ie(aF flv + 
fiFi iv ) XR dx tl dx v , with the minus sign due to the orienta- 
tion of the derivatives, and the contribution of the face 
at x R + dx x is ie(W~ 1 (aF ptu + [3F liU )W) {XR+dx x ) dx tl dx v , 



with Wi 



(x R +dx x ) 



1 



ieA x (xr) dx 



By Taylor 



expanding the second term, the joint contribution is 
ieD x (aF tlu + /3F ll ^) XR dx IJ 'dx u dx x , with no sums in the 
Lorentz indices. The contributions of the other two pairs 
of faces are similar, and the l.h.s. of (1) to lowest order 
is l + ie(D x [aF^ + fiF^] + cyclic perm.) XR dx fl dx u dx x . 
When evaluating the r.h.s. of (1) we can take the in- 
tegrand at any point of 51 since the differences are of 
higher order. In addition, the commutator term in J 
given in (4) is of higher order w.r.t. the first term in- 
volving the current. Therefore, the r.h.s. of (1) to lowest 
order is 1. + iej3 J ^ uX dx^ dx^ dx x . Equating the coefficients 
of a and /3 one gets the pair of (Hodge dual) Yang- Mills 
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equations. 

Let us discuss some consequences of (1). In order to 
write it for a given volume Q, we had to choose a ref- 
erence point xr on its border, and define a scanning of 
SI with surfaces and loops. If one changes the reference 
point and the scanning, both sides of (1) will change. 
However, the generalized non-abelian Stokes theorem (8) 
guarantees that the changes are such that both sides are 
still equal to each other. Therefore, one can say that (1) 
transforms "covariantly" under the change of scanning 
and reference point. In fact to be precise, the equation 
(1) is formulated not on SI but on the generalized loop 
space LVt = {7 : S 2 — > SI 1 north pole — > xr e dft}. The 
image of a given 7 is a closed surface E in Q contain- 
ing xr. A scanning of SI is a collection of surfaces E, 
parametrized by r, such that r = corresponds to the 
infinitesimal surface around xr and r = 2tt to <9S7. Such 
collection of surfaces is a path in LQ and each one cor- 
responds to SI itself. In order to perform each mapping 
7 we scan the corresponding surface E with closed loops 
starting and ending at xr, and each loop is parametrized 
by er, in the same way as we did in the arguments leading 
to (8). Therefore, the change of the scanning of SI corre- 
sponds to a change of path in LSI. In this sense, the r.h.s. 
of (1) is a path dependent quantity in LSI and its l.h.s. 
is evaluated at the end of the path. Of course, we do not 
want physical quantities to depend upon the choice of 
paths in LSI, neither on the reference point. Note that if 
we take, in the four dimensional space-time M, a closed 
tridimensional volume Sl c , then the integral Yang-Mills 
equation (1) implies that 



§ dCdrVJV- 



1 



(9) 



since the border 9S1 C vanishes, and the ordered integral 
of the l.h.s. of (1) becomes trivial. On the loop space 
LS7 C , S7 C corresponds to a closed path starting and end- 
ing at xr. Consider now a point 7 on that closed path, 
corresponding to a closed surface E, in such a way that 
Sli corresponds to the first part of the path and SI2 to the 
second, i.e. Sl c = Oi + O2 , and E is the common border of 
Sli and SI2 • By the ordering of the integration determined 
by (7) one observes that the relation (9) can be split 

f dQdrVJV- 1 f dCdrVJV- 1 

as P 3 e Jn 2 i"3e Jf! i = 1. However, by 

reverting the sense of integration along the path, one gets 
the inverse operator when integrating (7). Therefore, fli 
and S7JT 1 are two different paths (volumes) joining the 
same points, namely the infinitesimal surface around xr 
and the surface E, which correspond to their border. One 
then concludes that the operator P^e^n d i drV ^ v j s m _ 
dependent of the path, and so of the scanning of SI, as 
long as the end points, i.e. xr and the border 9S1, are 
kept fixed. 

The path independency of that operator can be used to 
construct conserved charges using the ideas of [2, 3] . First 



of all, let us assume that the space-time is of the form 
S x K, with 1R being time and S the spatial sub-manifold 
which we assume simply connected and without border. 
An example is when iS is the three dimensional sphere 

S 3 . It follows from (9) that Q s = P^s^^^ 1 = 1. 
That means that Qs is not only conserved in time, but 
also that there can be no net charge in S. In fact, there is 
the possibility of getting charge quantization conditions 
in such case (see [6, 7]). 

Let us now assume the space-time is not bounded, but 
still simply connected, like K 4 . We shall consider two 
paths (volumes) joining the same two points, namely the 
infinitesimal surface around xr, which we take to be at 
the time x° = 0, and the two-sphere at spatial infinity 
S'oo' 4 ', at x° = t. The first path is made of two parts. The 
first part corresponding to the whole space at x° = 0, i.e. 
the volume fl^} inside S&}°\ the two-sphere at spatial 
infinity at a; = 0. The second part is a hyper-cylinder 
S 2 ^ x /, where / is the time interval between a; = and 
x° = t, and S 2 ^ is a two-sphere at spatial infinity at the 
times on that interval. The second path is also made 
of two parts. The first one corresponds to the infinitesi- 
mal hyper-cylinder S$ x /, where S 2 is the infinitesimal 
two-sphere around xr and / as before. The second part 

corresponds to SlS, the whole space at time x° — t, i.e. 

2 (t) 

the volume inside Soo . From the path independency 
following from (9) one has that the integration of (7) 
along those two paths should give the same result, i.e. 
ViSl, x I) F(fi£ } ) = V(Cl£)V(S$ x I), where we have 
used the notation V (SI) = P 3 e^ d ( drV ^ v ^ anc [ wnere 
all integrations start at the reference point xr taken to 
be at x a — 0, and at the border S^ ^ of $7^ ■ In fact, 
one obtains V (SI) by integrating (7), and so one has to 
calculate JC = J ^ drV J V~ l , on the surfaces scanning 
the volume fl. We shall scan a hyper- cylinder S 2 x I 
with surfaces, based at xr, of the form given in figure 
(l.b), with t' denoting a time in the interval /. Each 
one of such surfaces are scanned with loops, labelled by 
r, in the following way. For < r < we scan the 
infinitesimal cylinder as shown in figure (l.a), then for 
X < T < x we scan the sphere S 2 as shown in figure 
(l.b), and finally for 4^ < r < 2ir we go back to xr with 
loops as shown in figure (l.c). The quantity K, can then 
be split into the contributions coming from each one of 
those surfaces as K = JC a + JCt, + IC C . In the case of the 
infinitesimal hyper-cylinder Sq x /, the sphere has in- 
finitesimal radius and so it does not really contribute to 
/Cft. We shall assume the currents and field strength van- 
ish at spatial infinity no slower than J M ~ 1/R 2+S , and 



F, 



l/m +d , with 6,5' > 0, for R 



Therefore 



the quantity J, given in (4), vanishes when calculated 
on loops at spatial infinity. Consequently, in the case of 
the hyper-cylinder S 2 ^ x /, the contribution to Kb coming 
from the sphere with infinite radius vanishes, and we have 
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that JC calculated on the surfaces scanning S 2 ^ x I and 
Sq x I is the same, and so V (S^ x /) = V (Sq x 7) . 
In fact there is more to it, since when we contract the 
radius of the cylinders in figure 1 to zero the loops in 
figures (l.a) and (l.c) become the same. Therefore, the 
quantities J calculated on them are the same except for 
a minus sign coming from the derivatives ^r, since the 
loops in figure (l.a) get longer with the increase of t, and 
in figure (l.c) the opposite occurs. In addition, the quan- 
tity V inside the the expression JC = dr V J V^ 1 
is insensitive to that sign since it is obtained by inte- 
grating (3) starting at xr in both cases. Therefore, 
it turns out that JC a + JC C = 0. The loops scanning 
the sphere in figure (l.b) have legs linking the reference 
point xr, at x° = 0, to the same space point but at 
x° = t' , i.e. x R . Therefore, when integrating (3) one gets 
Vx R = W{x R ,x R )- l V xt ^W{x R ,x R ), where W{x R ,x R ) is 
obtained by integrating (2) along the leg linking x R to 
x R , and where we have used the notation V x , meaning V 
obtained from (3) with reference point x. Using the same 
arguments and notation one obtains from (4) that, on the 
loops of figure (l.b), J XR = W(x R , xr)- 1 J xt > W(x R , x R ), 

and so JC b . XR = W(x R ,x R )~ 1 K, bx ^W(x t R ,x R ). The 

quantity V(ftoo) is obtained by integrating (7) and by 
scanning the volume with surfaces of the type 

shown in figure (l.b), and where the radius of S 2 
varies from zero to infinity keeping the point x R fixed. 
Therefore, from the above arguments one gets that 
= W{x t w x R )- 1 V xtR (n ( £)W(x t R ,x R ). One 
then concludes that such operator has an iso-spectral 
time evolution ^(0$) = U(t)V XR (*n$)U(t)- 1 , with 
U(t) = W(x R , x R )V (S$ xl). Therefore, its eigenval- 
ues, or equivalently Tr(V x< (tl£)) N , are constant in time. 
Note that from the Yang-Mills equations (1) one has that 
such operator can be written either as a volume or sur- 
face ordered integrals, and so we have proved (5). We 
have shown that, as a consequence of (9), such opera- 
tors are independent of the scanning of the volume. The 
reference point x R is on the border of the volume and 
so at spatial infinity. Then when we change the refer- 
ence point on the border to x R , the operator V x t (0^*) 
changes under conjugation by W(x R ,x R ). However, our 
boundary conditions implies that the field strength goes 
to zero at infinity and so the gauge potential is asymptot- 
ically flat, and consequently W(x R ,x R ) is independent 
of the choice of path joining the two reference points. 
Therefore, the conserved quantities are also independent 
of the base points. In addition, they are gauge invari- 
ant since, as shown in [5], V x t (fioo) — > gaYx* (^x)g R 1 , 
with g R being the group element, performing the gauge 
transformation, at x R . Note in addition that if V x t (f2^) 

has an iso-spectral evolution so does g c V x t (fiio ), with 
g c € Z{G). That fact has to do with the freedom we 



have to choose the integration constants of (3) and (7) 
to lie in Z(G), without spoiling the gauge covariance of 
(l)(see[5]). 




(a) (b) (c) 



FIG. 1: Surfaces of type (b) scan a hyper-cylinder S 2 x /. 

As an example consider a gauge theory for a gauge 
group G spontaneously broken to a subgroup H by 
a Higgs field 4> in the adjoint representation. For a 
BPS dyon solution one has Ei — sm9Di<fc, and Bi = 
COS0A0, with E i = -F M and B l = -\e ijk F jk , 
i,j,k = 1,2,3, with 9 being an arbitrary constant an- 
gle. At spatial infinity one has that D z <p — > j!^G(f), 
with f 2 = 1, r —> oo, and G{f) being an element 
of the Lie algebra of H, which is covariantly con- 
stant, i.e. D^G(r) — [8]. We have that the gauge 
field is asymptotically flat at spatial infinity, and so 
up to leading order one has = ^d^WW -1 , and 

so on 5^ (t) one has G(f) = WG R W'\ with G R 
being the value of G(f) at x R . Therefore, one has 

that V x t R {n { £) = P 2 e Js - ( ) L ^ J d " dT = 

exp [— ie{a cosO + (3 sm9)G R ]. Consequently, the con- 
served charges are given by the eigenvalues of G R , which 
contain among them the magnetic and electric charges of 
the dyon solution. Note that, even though we take G(f) 
at x R , the eigenvalues are independent of the choice of 
x R , since G(r) at different points at infinity are related 
by conjugation. 
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